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Abstract. Integrability, one of the classic issues in galactic dynamics and in general 
in celestial mechanics, is here revisited in a Riemannian geometric framework, where 
newtonian motions are seen as geodesies of suitable "mechanical" manifolds. The ex- 
istence of constants of motion that entail integrability is associated with the existence 
of Killing tensor fields on the mechanical manifolds. Such tensor fields correspond to 
hidden symmetries of non-Noetherian kind. Explicit expressions for Killing tensor 
fields are given for the N — 2 Toda model, and for a modified Henon-Heiles model, 
recovering the already known analytic expressions of the second conserved quantity 
besides energy for each model respectively. 
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1. Introduction 

The problem of integrability in classical mechanics has been very sem- 
inal. Motivated by celestial mechanics, it has stimulated a wealth of 
analytical methods and results. For example, the weaker request of 
only approximate integrability over finite times, or the existence of in- 
tegrable regions in the phase space of a globally non-integrable system, 
have led to the development of classical perturbation theory, with all 
its important achievements. However, deciding whether a given Hamil- 
tonian system is globally integrable or not still remains a difficult task, 
for which a general constructive framework is lacking. Besides its theo- 
retical interest, the problem of integrability is still relevant to a number 
of open problems among which we can mention a long standing one in 
galactic dynamics: the quest for the third integral of motion besides 
energy and angular momentum (Binney and Tremaine, 1987). In fact, 
the apparent absence of dynamical chaos in several models describing 
the motions of test stars in mean-field galactic gravitational potentials 
suggests that these models might be integrable. 
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The aim of the present paper is to draw attention to the rephrasing 
of this classical problem of integrability in a Ricmannian geometric 
language. Such a possibility exists because Hamiltonian flows can be 
identified with geodesic flows on Riemannian manifolds equipped with 
suitable metrics. The Riemannian geometric framework has hitherto 
proved very useful to tackle Hamiltonian chaos (Casetti, Pettini and 
Cohen, 2000), to understand the origin of chaos in non-Anosov flows 
(Pettini, 1993; Cerruti-Sola and Pettini, 1996; Cerruti-Sola, Pranzosi 
and Pettini, 1997; Pettini and Valdettaro, 1995; Casetti, Clementi and 
Pettini, 1996) and also to analytically compute Lyapunov exponents 
(Casetti, Clementi and Pettini, 1996). 

The existence of conservation laws, and of conserved quantities along 
the trajectories of a Hamiltonian system, is related with the existence 
of symmetries. The link is made by Nocthcr theorem (Arnold, 1978). 
A symmetry is seen as an invariance under the action of a group of 
transformations, and, in the case of continuous symmetries, this can be 
related also to the existence of special vector flelds, Killing vector fields 
on the mechanical manifold, generating the transformations. However, 
through Noetherian symmetries, and thus Killing vector fields, only a 
limited set of conservation laws can be accounted for. This is easily 
understood because only invariants that are linear functions of the 
momenta can be constructed by means of Killing vectors, while the 
energy, an invariant for any autonomous Hamiltonian system, is already 
a quadratic function of the momenta. The possibility of constructing 
invariants along a geodesic flow, that are of higher order than linear 
in the momenta, is related with the existence of Killing tensor fields 
on the mechanical manifolds (Eisenhart, 1964). In the present paper 
we discuss all these facts and we show how it is possible to explicitly 
work out the components of the Killing tensors associated with two 
integrable models: an A?^ = 2 Toda model and a modified version of the 
Hcnon-Heiles model, and hence to obtain the analytic expressions of 
the second integral of motion besides energy. 

In general, the components of any Killing tensor field on a me- 
chanical manifold are solutions of a linear, non-homogeneous system 
of first order partial differential equations. As the number of these 
equations always exceeds the number of the unknowns, the system is 
always overdetermined. The existence of Killing tensors thus requires 
compatibility. However, compatibility is generically very unusual, hence 
a possible explanation, at least of qualitative kind, of the exceptionality 
of integrability with respect to non-integrability. 

For the sake of clarity and self-containedness, we briefiy recall some 
basic points about the geometrization of newtonian mechanics, about 
Killing vector fields and Killing tensor fields in Sections II, HI, and IV, 
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respectively. Section V contains the original results mentioned above 
about the relationship between integrability and the existence of Killing 
tensor fields on the mechanical manifolds. 



2. Geometric formulation of Hamiltonian dynamics 

Let us briefly recall a few basic points about the geometrization of 
newtonian dynamics in a Riemannian geometric framework. This ap- 
plies to dynamical systems described by standard Hamiltonians, i.e. 
H{p, q) = ^a^-' {q)piPj + V (q) , with the shorthands p = (pi, . . . ,pn) and 
q = {qi,--- : 9iv)- Equivalently, we can describe these systems through 
Lagrangian functions L{q,q) = ^aij{q)q^q^ — V{q). 

According to Maupertuis' principle of stationary action, among all 
the possible isoenergetic paths 7(t) with fixed end points, the paths 
that make vanish the first variation of the action functional 




are natural motions. 

The kinetic energy is a homogeneous function of degree two, 
hence 2W = qidL/dqi , and Maupertuis' principle reads 

SA = S [ 2Wdt = . (2) 

hit) 

The configuration space M of a system with degrees of freedom is an 
A-dimensional differentiable manifold and the lagrangian coordinates 
(gi, . . . , qisf) can be used as local coordinates on M. The manifold M is 
naturally given a proper Riemannian structure. In fact, by introducing 
the matrix 

gij = 2[E - V{q)]aij (3) 

Eq.(2) becomes 

6 f 2Wdt = S [ (gnq'q^Y^^ dt = 6 f ds = , (4) 

J-/{t) J fit) ^ ^ Jfis) 

so that the newtonian motions fulfil the geodesic condition on the man- 
ifold M, provided we define ds as its arclength. The metric tensor gj of 
M is defined through its components by Eq.(3). This is known as Jacobi 
(or kinetic energy) metric. Denoting by V the canonical Levi-Civita 
connection on (M, gj ) , the geodesic equation 

V^7 = (5) 
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becomes, in the local coordinates {q^, . . . , q^), 

^+n,^^ = 0, (6) 
ds^ ds ds ^ ^ 

where the Christoffel coefEcients F*-^ are the components of V defined 
by ^ 

Tjfe = ^g^"^ {dj9km + dk9mj - dmgjk) (7) 

= -^[SidjV + 5idkV - diVa^'ajk] + ^a*^[5«a,fc + dkaij - diajk] , 

where di = d/dq^. Without loss of generality consider gij = 2[E — 

V{q)\5ij, so that 

^ik = -^mv + SpkV - Sjkd'V] . (8) 
and from Eq. (6) we get 

d'^q' 1 



ds2 2{E - V) 



d{E-V)dq^ dq^ ,^d{E-V) dq^ dq"'' 
2 Z 1 — 1 9 5 9ki 



dqj ds ds dqj ds ds 

and, using ds"^ = 2{E — V)"^ dt"^, these equations finally yield 
d^q' dV 



= 0, 
(9) 



dt'^ dqi 
which are Newton equations. 



z = l,...,iV. (10) 



3. Killing vector fields 

On a Riemannian manifold, for any pair of vectors V and W , the 
following relation holds 

where (y,W) = gijV'W^ and V/ds is the covariant derivative along a 
curve 7(s). If the curve 7(5) is a geodesic, for a generic vector X we 
have 

^(X,7) = (^.7> + (X.^) = (^.7) ^ (V,X,,) (12) 



where (V-yX)' = ^^|r + ^]k^-^''' ^° components it reads 

^^{XiV^) = v^Vi{XjV^) (13) 
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where = dx^ /ds; with Xjv'^ViV^ = XjVjj^ = - because geodesies 
are autoparallel - this can be obviously rewritten as 



(Xiv') = -v^v\ViXj + VjXi) (14) 



telhng that the vanishing of the l.h.s., i.e. the conservation of XjV* along 
a geodesic, is guaranteed by the vanishing of the r.h.s., i.e. 

V(iXj-) = ViXj + VjXi = , i,j = 1,..., dimME. (15) 

If such a field exists on a manifold, it is called a Killing vector field 
(KVF). Equation (15) is equivalent to Cxd = 0, where C is the Lie 
derivative. On the mechanical manifolds {Me, gj), being the unit vector 
^ - tangent to a geodesic - proportional to the canonical momentum 
Pk = = q'^, {cLij = %)) the existence of a KVF X implies that the 
quantity, linear in the momenta, 

(16) 

is a constant of motion along the geodesic flow. Thus, for an N degrees 
of freedom Hamiltonian system, a physical conservation law, involving 
a conserved quantity linear in the canonical momenta, can always be 
related with a symmetry on the manifold {ME,gj) due to the action 
of a KVF on the manifold. These are conservation laws of Noetherian 
kind. The equation (15) is equivalent to the vanishing of the Poisson 
brackets 

the standard definition of a constant of motion J{q,p). In fact, a linear 
function of the momenta 

J{q,p) =YlCi{q)Pi , (18) 

i 

if conserved, can be associated with the vector of components 

Xk = [E-V{q)]Ck{q). (19) 

The explicit expression of the system of equations (15) is obtained by 

writing in components the covariant derivatives associated with the 
connection coefficients (8) and it finally reads 



[E-V{q)] 



dCiiq) ^ dC,{q) 



dq^ dq^ 



N 



8V 

E a^^^(^) = 0' (20) 
fe=i " 
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or equivalently 



which, according to the principle of polynomial identity, yields the 
following conditions for the coefficients Ci{q) 



dq' 

N 



1,...,N (22) 



One can easily check that the same conditions stem from Eq.(17). 
As an elementary example, we can give the explicit expression of the 
components of the Killing vector field associated with the conservation 
of the total momentum P{q,p) = J2k=iPk ■ 

In this case the coefficients are Ci{q) = 1, so that the momentum 
conservation can be geometrically related with the action of the vector 
field of components 

Xi = E-V{q), i = l,...,N (23) 

on the mechanical manifold. At least this class of invariants has a 
geometric counterpart in a symmetry of {ME,gj). 

However, in order to achieve a fully geometric rephrasing of intc- 
grability, we need something similar for any constant of motion. If 
a one-to-one correspondence is to exist between conserved physical 
quantities along a Hamiltonian flow and suitable symmetries of the 
mechanical manifolds {ME,gj), then integrability will be equivalent to 
the existence of a number of symmetries at least equal to the number 
of degrees of freedom (= dim Me)- 

If a Lie group G acts on the phase space manifold through completely 
canonical transformations, and there exists an associated momentum 
m,apping^ , then every Hamiltonian having G as a symmetry group, with 
respect to its action, admits the momentum mapping as constant of 
motion (Abraham and Marsden, 1987). These symmetries are usually 



^ This happens whenever this action corresponds to the Ufting to the phase space 
of the action of a Lie group on the configuration space. 
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referred to as hidden symmetries because, even though their existence 
is ensured by integrabihty, they are not easily recognizable^ . 




4. Killing tensor fields 

Let us now extend what has been presented in the previous section 
about KVFs, trying to generalize the form of the conserved quantity 
along a geodesic flow from J = XiV^ to J = Kjj^j2...jj.v^^v^^ . . . v^'' , with 
^h32--jr * tensor of rank r. Thus, we look for the conditions that entail 

-^^{Kj,j,,„jy^v^^ . . . = v^Vj{Kj,j,,„jy^v^^ . . . v^^) = . (24) 

In order to work out from this equation a condition for the existence of 
a suitable tensor ETj^jj-.-ir) which is called a Killing tensor field (KTF), 
let us first consider the 2r rank tensor Kj^j^^^^j^v'^'^v'^^ . . .v'^^ and its 
covariant derivative along a geodesic, i.e. 

= v'H'^...v'r^^VjKj,j,,„j^ (25) 

where we have again used v^VjV^'' = along a geodesic, and a standard 
covariant differentiation formula (Doubrovine et al., 1979). Now, by 
contraction on the indices ik and jk the 2r-rank tensor of the r.h.s. of 
Eq.(25) provides a new expression for the r.h.s. of Eq.(24) which reads 

^^{Knn...jy'v^' . . . = vV . . . v^^v^V ^K^,^,...^^) , (26) 

where V ^jKj^j^,,,^^) = V,i^ju2...jV + ^.nKjj^...jr + ■■■ + ^jrKnn...j,_,j, 
as it can be easily understood by rearranging the indices of the summa- 
tions in the contraction of the 2r-rank tensor in the last part of Eq. (25); 

^ An interesting account of these hidden symmetries can be found in (Olshanetsky 

and Pcrclomov, 1981) where it is surmised that integrable motions of A'^ degrees of 
freedom systems are the "shadows" of free motions in symmetric spaces (for example 
euclidean spaces M", hyperspheres S", hyperbolic spaces H") of sufficiently large 
dimension n> N. 
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(a direct check for the case N = r = 2 is immediate). The vanishing 
of Eq.(26), entailing the conservation of ^jij2-.-ir^"'^^''^ • • • '""''^ along a 
geodesic flow, is therefore guaranteed by the existence of a tensor field 
fulfilling the conditions 

Vo-^,U2...>)=0' (27) 

these equations generalize Eq.(15) and give the definition of a KTF on 
aRiemannian manifold. These N'^'^^ equations in {N+r — l)\/r\{N—l)\ 
unknown independent components^ of the Killing tensor constitute an 
overdetermined system of equations. Thus, a-priori, we can expect that 
the existence of KTFs has to be rather exceptional. 

If a KTF exists on a Riemannian manifold, then the scalar 



is a constant of motion for the geodesic flow on the same manifold. 

Let us consider, as a generalization of the special case of rank one 
given by Eq.(18), the following constant of motion 

Jiq,p)= J2 C'm2...rfP2 •••f'i^, (29) 
{ii,i2,...,jjv} 

which, with the constraint ii + ^2 + • • • + = r, is a homogeneous 
polynomial of degree r. The index ij denotes the power with which the 
momentum pj contributes. If r < then necessarily some indices ij 
must vanish. By repeating the procedure developed in the case r = 1, 
and by identifying 

J. . ^ dq^^dq^^ dq^^ 

J{q,p)=Kj,j,..,^^^...^ (30) 

we get the relationship between the components of the Killing tensor 
of rank r and the coefficients Ci^i^^^ij^ of the invariant J{q,p), that is 

--2'/^[E-Viq)Ya,,,...,^ . (31) 




With the only difference of a more tedious combinatorics, also in this 
case it turns out that the equations (27) are equivalent to the vanishing 
of the Poisson brackets of J{q, p) , that is 

{^,-^} = O^V(,X,,,,...,-^) = 0. (32) 

^ This number of independent components, i.e. the binomial coefficient (^^^~^), 
is due to the totally symmetric character of Killing tensors. 
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Thus, the existence of Kilhng tensor fields, obeying Eq.(27), on a me- 
chanical manifold {M,gj) provide the rephrasing of intcgr ability of 
Newtonian equations of motion or, equivalcntly, of standard Hamilto- 
nian systems, within the Riemannian geometric framework . 

At first sight, it might appear too restrictive that prime integrals 
of motion have to be homogeneous functions of the components of p. 
However, as we shall discuss in the next Section, the integrals of motion 
of the known integrable systems can be actually cast in this form. This 
is in particular the case of total energy, a quantity conserved by any 
autonomous Hamiltonian system. 



5. Explicit KTFs of known integrable systems 



The first natural question to address concerns the existence of a KT 
field, on any mechanical manifold {M,gj), to be associated with total 
energy conservation. Such a KT field actually exists and coincides with 
the metric tensor gj, in fact it satisfies^ by definition Eq.(27). 

One of the simplest case of integrable system is represented by a 
decoupled system described by a generic Hamiltonian 



N 



N 



(33) 



for which all the energies Ei of the subsystems Hi, i = I, . . . ,N, are 
conserved. On the associated mechanical manifold, N KT fields of rank 
2 exist, they are given by 

Kfi = 6,kmq,)[E - V{q)] + 6}[E - Viq)]'} . (34) 

In fact, these tensor fields fulfil Eq.(27) which explicitly reads 



(i) 
Im 



dK 



mk 



dK, 



(i) 
kl 



(35) 



dq'^ ' ' dq'^ kl-^'-jm km jl ~ Im jk 

k,l,m = I, . . . , N . 

The conserved quantities J^^\q,p) are then obtained by saturation of 
the tensors with the velocities dq/ds 

jk = l dS dS ^ r y^j ^ 

(36) 

* A property of the canonical Levi-Civita connection, on which the covariant 
derivative is based, is just the vanishing of S7g. 



EPk ,Pi _ p 
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This equation suggests that to require that the constants of motion have 
to be homogeneous polynomials of the momenta is not so restrictive as 
it might appear, in fact, through the following constant quantity 

1 N 2 



homogeneous of second degree in the momenta, any even degree poly- 
nomial of the momenta can be made homogeneous. The possibility of 
inferring the existence of a conservation law from the existence of a 
KTF on (M,gj) is thus extended to the constants of motion given by 
a sum of homogeneous polynomials whose degrees differ by an even 
integer 

J(p, q) = p(0 (p) + p('--2) (p) + . . . + (38) 
+p(r-^-)(p)eC'^{q)\pi,...,PN] 

T 

homdegP^ = s s = r,r — 2, . . . ,r — 2[-] 
so that it can be recast in the homogeneous form 

1 N 2 



+p('-2")(p) 



1 N 2' 



5.1. NONTRIVIAL INTEGRABLE MODELS 

Nontrivial examples of nonlinear integrable Hamiltonian systems are 
provided by the following Hamiltonians 

N ( 2 \ 



i=l 



known as the Toda model (Toda, 1970), which is integrable for any 
given pair of the constants a and 6, and 

^ = E| + ^ E^n -EM^ (41) 

i=l \i=\ J i=\ 

which is completely integrable for any Ai, . . . , Aat (Choodnovsky and 
Choodnovsky, 1978). Recursive formulae are available for all the con- 
stants of motion of the Toda model at any N (Henon, 1974), and 
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also for the second Hamiltonian the exact form of first integrals is 
known (Choodnovsky and Choodnovsky, 1978). In both cases, the first 
integrals are polynomials of given parity of the momenta so that, on 
the basis of what we have said above, each invariant J*^*\ i = 1, . . . , N 
can be derived from a KTF on {M,gj). Thus, integrability of these 
systems admits a Riemannian-geometric interpretation. 



5.2. The special case of N=2 Toda model 



Let us consider the special case of a two-degrees of freedom Toda model 
described by the integrable Hamiltonian^ 

H = \{pI +pI) + Ya [e''^'"^" + e''"'^" + e"""] - \ ■ (42) 

Prom what is already reported in the literature (Henon, 1974), we 
know that a third order polynomial of the momenta has to be found 
eventually, therefore we look for a rank-3 KT fulfilling 

ViKjki + VjKiki + VkKiji + ViKijk = , i,j,k,l = 1,2 (43) 

where, associating the label 1 to x and the label 2 to y, {{i,j,k,l)} = 
{(1, 1, 1, 1); (1, 1, 1, 2); (1, 1, 2, 2); (1, 2, 2, 2); (2, 2, 2, 2)}. The computation 
of the Christoffel coefficients according to Eq.(8) yields 

1 _ -d^V 1 _ dccV 2 _ 9yV 

ill — _ r _ TTT rr 1 J- oo — _ r _ TTZ rr •, i i 



2[E-V{x,y)]' 2[E-Vix,y)]' 2[E-Vix,y)]' 

J- 99 — _ r _ TTT TT , J- 19 — _ r _ TTT 7T , i 



22 2[E-V{x,y)]' 12 2[E-Vix,y)]' 2[E-Vix,y)]- 

(44) 

Prom Eq.(43) we get the system 

Vii^iii = 

V1K122 + V2i^ll2 = 

V2KU1 + 3V1K211 = 

V1K222 + 3V2i^l22 = 

V2i^222 = (45) 

whence 

d^cKui - Sriii^in - 3rj^K2n = 

^ This is derived from an A'^ = 3 Hamiltonian (40) by means of two canoni- 
cal transformations of variables removing translational invariance, see for example 
(Boccaletti and Pucacco, 1996); the third order expansion of this new Hamiltonian 
yields the Henon-Heiles model of Eq.(49) with C = D = 1. 
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4rii^212 — ^22Kui ~ r|2^211 = 

dyKui + Sd^K2n - GTl^Kiu - erl2Kn2- 

Gr\^K2ii - 6r?iir2i2 = 

d^K222 + SdyKi22 - 6rliKi22 " Gr^^ 7^222" 

6r^2^112 ~ 6r22i^212 = 
dyK222 - 3r^2^122 - 3r^2^222 = (46) 

with the ChristofFcl coefficients given by Eq.(44), where one has to 
replace V{x,y) with the potential part of the Hamiltonian (42) and 
dxV, dyV with its derivatives. The general method of solving a linear, 
non-homogeneous system of first-order partial differential equations in 
more than one dependent variables is sketched in Appendix. However, 
finding the explicit solution to the system of equations (46) is much 
facilitated because we already know a-priori that this system is com- 
patible and thus admits a solution, and we also have strong hints about 
the solution itself because the general form of the integrals of the Toda 
model is known (Henon, 1974). The KTF, besides the metric tensor, 
for the model (42) is eventually found to have the components 

Km = 2{E-Vf[3dyV + 4{E-V)] 

= 8{E - Vf + 1{E- l/)2[e2j/-2V3x + ^2y+2V3x __ 2g-4y] 

i^i22 = 2{E - Vf[dyV - AiE - V)] 

= -2A{E -Vf + ^{E- y)2[e22/-2V3x ^ ^2y+2V3x _ c^^-iy^ 

Ku2 = -2{E - Vfd^V = ^{E- y)2(e2i/+2^^ _ e2y-2V3x) 

6 

K222 = -^{E - Vfd^V = ^{E- y)2(e2i/+2^- _ ^2^-2^3.)^ (47) 

as can be easily checked by substituting them into Eqs.(46). Hence, the 
second constant of motion, besides energy, is given by 

, \ _ ^1'^ _ tr ^1'' ^ 

"J y^^y^PxiPy) — ijk~l '1 3 ~ ^ijk' 



ds ds ds ^ dt dt dt 2^/2[E - V{x,y)]^ 

" 2^/2[E-V(a;,y)]3 ^^'''^' ^ 3ifi22PxP^ + SKu2plPy + ^^222^^) 

= 8p,{pl - 3pI) + (p, + V3py)e^y-^-^'- - 2p,e-^y + (p, - V3py)e'y+^^'^ 

(48) 
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which coincides with the expression already reported in the hterature 
(Lieberman and Lichtenberg, 1992) for the Hamiltonian (42). 

5.3. The generalized Henon-Heiles model 

Let us now consider the two-degrees of freedom system described by 
the Hamiltonian 

H = 1{pI+pI) + 1 {x' + y') + Dx'y -Icy'. (49) 

This model, originally derived to describe the motion of a test star in an 
axisymmetric galactic mean gravitational field, provided one of the first 

numerical evidences of the chaotic transition in nonlinear Hamiltonian 
systems (Henon and Heiles, 1964). Henon and Heiles considered the 
case C = D = 1. The existence of a chaotic layer in the phase space 
of this model means lack of global intcgrability. However, by means 
of the Painlcvc method, it has been shown (Chang, Tabor and Weiss, 
1982) that for special choices of the parameters C and D this system is 
globally integrable. Let us now tackle intcgrability of this model from 
the viewpoint of the existence of KT fields on the manifold {M,gj). 
We first begin with the equations for a Killing vector field. By means 
of Eqs.(22) we look for possible coefficients Ci{x,y), C2{x,y), thus 
obtaining 

C, = Ci(9), C2 = C2(x) 

^ + !?^ = (50) 

dy dx 

x(l + 2Dy)Ci{y) + {y + Dx^ - Cy^)C2{x) = 
From the second equation of Eqs.(50) it follows that 

dCi{y) dC2{x 



dy dx 



= cost. (51) 



whence, denoting with a the constant, the possible expressions for 
Ci{y) and C2{x) are only of the form Ci{y) = —ay +(5 , C2{x) = ax+7, 
that, after substitution into the last equation of Eqs. (50), imply 

{x + 2Dxy){-ay + /?) + (y + Dx^ - Cy^){ax + 7) = 0, (52) 

which has only a non-trivial solution for C = D = Q. On the other 
hand, for these values of the parameters the potential simplifies to 
V{x^ y) = + whence the existence of the Killing vector field X 
of components Xi = y and X2 = —x which is due to the invariance 
under rotations in the xy plane. 
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Let US now consider the case of a rank-2 KTF. Equations (43) 
become 

ViKjk + VjKik + VkKij = , i,j,k = 1,2 (53) 

where, associating again the label 1 to x and the label 2 to y, {(i, j, k)} = 
{(1, 1, 1); (1, 1, 2); (1, 2, 2); (2, 2, 2)}. The ChristofFel coefficients are still 
given by Eq. (44) , where we have to use the potential part of Hamilto- 
nian (49) so that dxV{x,y) = x+2Dxy and dyV{x,y) = y+Dx'^ — Cy'^. 
The KTF equations are then 

ViKii = 
2V1K12 + V2K11 = 
V1K22 + 2V2K12 = 

V2K22 = (54) 

whence 

25,Ki2 + dyKn 

d^K22 + 2dyKi2 

(55) 

Since the Hamiltonian (49) is not integrable for a generic choice of 
the parameters C and D, we can reasonably expect that the generic 
property of the above overdetermined system of equations is incompat- 
ibility, i.e. only the trivial solution Kij = exists for the overwhelming 
majority of the pairs {C,D). However, the existence of special choices 
of C and D for which the Hamiltonian is integrable suggests that this 
overdetermined system can be compatible in special cases. For example, 
when D = the variables x and y in (49) are decoupled and thus two 
KT fields of rank 2 exist according to Eq.(34). 

A non trivial solution for the system (55) must exist at least for the 
pair (C = — 6, D = 1). In fact, in this case the modified Henon-Heiles 
model is known to be integrable (Chang, Tabor and Weiss, 1982). An 
explicit solution for the system (55) is eventually found to be given by 

Kii = (3 - 4.y){E - V{x, y)f + x\x' + Ay"" + 4y + 3)(E - V{x, y)) 

Ki2 = 2x{E-V{x,y)) (56) 

K22 = \i^^ + + % + 3) (S - V{x, y)) . 



d^Kn - 2riiKn - 2r?iK2i = 

- iT\^Kn - i^Tl^ + 2T\^)Ku - 2TI^K22 = 

-2Tl^K^^-{AT\^ + 2Tl^)K^2-ATl^K22 = 

dyK22 - 2r^2^i2 - 2ri2i^22 = 
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The associated constant of motion is therefore 



J{x,y,Px,Py) = 



1 



(KupI + 2KuPxPy + K22pI) = 



{E-V{x,y)f 



= + 4x^2/^ - ply + 4:PxPyX + 4x^y + 3pl + . 



(57) 



This expression is identical to that reported in (Chang, Tabor and 
Weiss, 1982), worked out for the same values of C and D with a 
completely different method based on the Painleve property. 



Let us now summarize the meaning of the results presented above and 
point out the open problems. 

— Besides qualitative and quantitative descriptions of chaos, within 
the framework of Riemannian geometrization of newtonian me- 
chanics also integrability has its own place. The idea of associating 
KTFs with integrability is not new, though this has been essen- 
tially developed in the context of classical General Relativity, see 
for example (Baleanu, 1998; Gibbons et al., 1993; Sommers, 1973) 
and references quoted therein. Recently, also an extension to clas- 
sical newtonian mechanics has been considered in (Rosquist and 
Pucacco, 1995), where integrability conditions for quadratic invari- 
ants were obtained, and where the authors concluded saying "It is 
of considerable interest to develop these techniques and use them 
to look for fixed energy invariants of physically interesting models 
such as the Henon-Heiles potential and others": this is just what 
has been done in our present work. 

— The reduction of the problem of integrability of a given Hamilto- 
nian system to the existence of suitable KTFs on (MejQj) offers 
several reasons of interest, in particular we have seen that the 
system of equations in the unknown components of a KTF of a 
preassigned rank is overdetermined, thus - at a qualitative level - 
integrability seems a rather exceptional property, and the larger N 
the "more exceptional" it seems to be, because of the fastly growing 
mismatch between the number of unknowns and the number of 
equations. In principle, the existence of compatibility conditions 
for systems of linear, first-order partial differential equations could 
allow to decide about integrability prior to any explicit attempt at 
solving the equations for the components of a KTF. Even better. 



6. Concluding remarks 
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there are geometric constraints to the existence of KTFs, early 
results in this sense are reported in (Yano and Bochner, 1953), 
so that it seems possible, at least in some cases, to devise purely 
geometric criteria of non-integrability. For example, hyperbolicity 
of compact manifolds excludes (Yano and Bochner, 1953) the exis- 
tence of KTFs, and this is consistent with the property of geodesic 
flows on compact hyperbolic manifolds of being strongly chaotic 
(Anosov flows). 

— In the present paper, before working out the second invariant 
besides energy for two integrable models, we already knew that 
a KTF had to exist and of which rank (because of the degree 
of the polynomial invariant), thus we knew that the system of 
equations to be solved was compatible. Whereas, in general we 
lack a criterion to restrict the search for KTFs to a small interval 
of ranks, and this constitutes a practical difficulty. Nevertheless, 
since the involution of two invariants translates into the vanishing 
of special brackets - the Schouten brackets (Sommers, 1973) - 
between the corresponding Killing tensors, a shortcut to prove 
integrability, for a large class of systems fulfilling the conditions of 
the Poincare- Fermi theorem (Poincare, 1892; Fermi, 1923), might 
be to find only one KTF of vanishing Schouten brackets with the 
metric tensor. In fact, for quasi-integrable systems with > 3, 
the Poincarc-Fermi theorem states that generically only energy is 
conserved, thus if another constant of motion is known to exist 
(apart from Noetherian ones, like angular momentum) then the 
system must be integrable and in fact there must be N constants 
of motion. 

— At variance with Killing vectors, which are associated with Noethe- 
rian symmetries and conservation laws, Killing tensors no longer 
have a simple geometrical interpretation (Gibbons et al., 1993; 
Rosquist, 1989), therefore the associated symmetries are non-Noetherian 
and hidden. 

The present paper contributes the subject of a Riemannian geometric 
approach to integrability with constructive examples that non-trivial 
constants of motion besides energy can be derived from KTFs for two 
degrees of freedom integrable Hamiltonians of physical interest. This 
approach to integrability deserves further attention and investigation. 
In fact, among the other reasons of interest, by considering, for ex- 
ample, the standard Henon-Heiles model (C = £) = 1), we could 
wonder whether the regular regions of phase space correspond to a 
local fulfilment of the compatibility conditions of the system (55) , this 
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would lead to a better understanding of the relationship between geom- 
etry and stability of newtonian mechanics. Moreover, we could imagine 
that, by suitably defining weak and strong violations of these compat- 
ibility conditions, we could better understand the geometric origin of 
weak and strong chaos in Hamiltonian dynamics (Pettini and Landolfi, 
1990; Pettini and Ccrruti-Sola, 1991) and, perhaps, this might even 
suggest a starting point to develop a "geometric perturbation theory" 
complementary to the more standard canonical perturbation theory. 

Finally, the celebrated problem of the third integral in galactic dy- 
namics could find here a new constructive, and hopefully useful, ap- 
proach. 
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Appendix 

Let us now briefly sketch a classical method (Forsyth, 1959) of solving 
systems of linear, first-order, partial differential equations in several 
dependent variables, denoted hy zi, . . . , Zm, and two independent vari- 
ables, denoted by x and y. Writing Xi = dzi/dx and Yi = dzi/dy, 
equations (46) and (55) are in the form 

m 

Xi = Bi + J2Ai,Ys (58) 

s=l 

with an obvious meaning of the coefficients Ais and S^; i = 1, . . . ,m. 
The first step consists of "diagonalizing" the above system, writing 
equivalent first order equations, or systems of equations, in only one 
dependent variable. Thus equations (58) are multiplied by Ai, . . . , A^, 
a set of multipliers, and summed to give 

m m m 

^ XiBi - XiXi + J2 ^i^isYs = . (59) 

1=1 i=l j,s=l 

Since the derivatives must fulfil the obvious relations 

dzi - Xidx - Yidy = , (60) 
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by combining Eqs.(59) and (60) the following system of ordinary equa- 
tions is formed 

s = 1, . . . ,m, whence, by putting dy = jidx , we obtain 

m 

J2 XiAs + I^Xs = (62) 
1=1 

so that, by solving the critical equation det(^ + /jI) = 0, where I is 
the identity matrix, the numbers Ai, . . . , Am can be eliminated among 
these equations obtaining a set of ratios Ai : A2 : • • • : A^, then the 
quantities defined by 

= o (63) 

are uniquely determined. The above system of ordinary equations is 
finally rewritten as 

aidzi + . . . + amdzm = dx 

dy = fidx (64) 

i.e. two linear characteristic equations for each root of the criti- 
cal equation. Now, if u{x,y, zi, . . . , Zm) = const is an integral of the 
equations (64), then it fulfils also the system of m equations 

du du du 

\- ai— + ain— = 

ozi ox oy 



du du du ^ 

+ am^ + am/^^ = (65) 



ozm ox dy 

since i = 1, . . . ,m, there are m equations in m -|- 2 variables; an in- 
tegral of this system, involving any of the dependent variables, is an 
integral of the original system of equations, and if each of the roots of 
the critical equation leads to an integral, then the ensemble of these 
integrals provides an integral equivalent of the original system. With 
the substitutions 

Zi = Wi 

Zi = Ci + {wi- ci)wi , i = 2,...,m (66) 
the following single equation is constructed 

du du du , , 

Ui^ + U2^ = (67) 



dwi dx dy 
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where Ui = ai +J2iL2 cui'^i a-^d U2 = Oiiii+Y^^2 l^onwi. Finally, the in- 
tegration of the equation (67) proceeds by integrating the characteristic 
equations 

= - = j| (68) 

taking W2,---, as non- varying quantities, and then proceeding in a 
standard way (Courant and Hilbert, 1962). 
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